Experimental evidence for the intricate free-energy landscape for a soft glassy system 
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In the free energy landscape picture of glassy systems, the slow dynamics characteristic of these 
systems is believed to be due to the existence of a complicated free-energy landscape with many 
local minima. We show here that for a colloidal glassy material multiple paths can be taken through 
the free energy landscape, that can even lead to different 'final' non-ergodic states at the late stages 
of aging. We provide clear experimental evidence for the distinction of gel and glassy states in the 
system and show that for a range of colloid concentrations, the transition to non-ergodicity can 
occur in either direction (gel or glass), and may be accompanied by 'hesitations' between the two 
directions. This shows that colloidal gels and glasses are merely global free-energy minima in the 
same free energy landscape, and that the paths leading to these minima can indeed be complicated. 

PACS numbers: 



The main issue in the study of glasses is that from the 
point of view of their molecular structure, they closely 
resemble liquids. On the other hand, their mechani- 
cal properties are much closer to those of solids: ordi- 
nary window glass, for instance, does not flow on human 
timescales To explain the extremely high viscosity 
of glassy systems, it is commonly accepted that the mo- 
tion of molecules or particles that constitute the glass 
are blocked by the neighboring molecules, who in turn 
are blocked also by their neighbors and so on, making it 
impossible for the system to flow. 

Translated in terms of the free energy of the system, 
the paradigm for glasses is that of a complicated free en- 
ergy landscape, with many local minima of the free en- 
ergy [2, 3]. The blocking of molecules or particles is then 
equivalent to saying that they are trapped in a local min- 
imum of the free energy from which they cannot escape. 
This is most evident in the " aging" of such systems 0, ^ : 
just after a quench into the glassy phase, the system still 
evolves in the sense that both the mechanical properties 
and the diffusion coefficient change in time. The inter- 
pretation of this time evolution is that at early times 
after its formation, the system is able to access at least 
part of the phase space, and can get out of local minima 
by thermal activation. However, as time goes on, the 
system finds deeper and deeper minima, more difficult 
to escape from, and consequently the evolution becomes 
slower. Because of this, the system cannot reach ther- 
modynamic equilibrium: it remains non-ergodic. During 
this aging, the viscosity increases and the diffusion co- 
efficient of the particles decreases, emphasizing the link 
between the blocking of the particle motion and absence 
of flow. 

However, although this provides an appealing intuitive 
picture of glassy dynamics (or rather, the absence of dy- 
namics), to our knowledge there is no direct experimen- 
tal evidence for the existence of such a complicated free- 
energy landscape with many local minima [3]. In this 



Letter we are able to provide such evidence by showing 
that for a soft glassy material multiple paths can be taken 
through the free energy landscape, that can even lead to 
different 'final' non-ergodic states at the late stages of 
aging. 

The system we consider is a suspension of anisotropic 
and charged colloidal particles suspended in water: 
Laponite, a synthetic clay 0, @, S B H, K [HI- The 
study of colloids allowed for a significant contribution to 
elucidating the basic physics of glass transition [l^, [l^ ■ 
In colloidal systems, as the particle volume fraction is in- 
creased, the particles become increasingly slower and for 
even higher volume fractions the glass transition is en- 
countered. On the other hand, colloidal gels are known 
to form at extremely low volume fractions 10^^ — 10^^ 
in the presence of strong attractions [3]. Gelation and 
the glass transition have important similarities. Both 
are ergodic to non-ergodic transitions that are kinetic, 
rather than thermodynamic in origin, and distinguishing 
between these two types of non-ergodic states experimen- 
tally is a longstanding controversy [lit . The experiments 
reported below provide direct criteria for distinguishing 
gels from glasses. This allows us to show that for a range 
of Laponite concentrations, two distinctly different non- 
ergodic states can result at late times: either the glass 
or the gel forms at late times with roughly equal prob- 
ability. There is no way to tell beforehand which of the 
two options will be taken by the sample, suggesting that 
there are at least two metastable minima in the system. 
In addition, our results show that the free energy land- 
scape is indeed complicated, since a number of samples 
are observed to hesitate between the two options for a 
long time and an initial evolution in one of the two di- 
rections can lead to a final state that is the other one. 

The formation of non-ergodic states in our sys- 
tem is followed in time using light scattering. The 
non-crgodicity parameter (Debye- Waller factor) deduced 
from these measurements, (which quantifies the fraction 
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of frozen-in fluctuations pj|) falls onto either of two mas- 
ter curves, showing that two possible routes towards non- 
ergodicity exist. Measurement of translational diffusion 
and scattered intensity then allows to unambiguously 
identify these two distinct non-ergodic states as colloidal 
gel and glass. 

For the experiments, we prepare Laponite XLG sus- 
pensions in ultrapure water. After stirring for two hours, 
we filtered using Millipore O.S^m filter units to obtain 
a reproducible initial state This instant defines the 
zero of aging time, ta = 0. A standard dynamic light 
scattering (A = 632. Snm) measures the time-averaged in- 
tensity correlation functions gtiq,t) — ^''^S'^r/'''^nw^^' i 
scattering wave vector q = ^^sin(|), in which 6 = 90° 
is the scattering angle. The correlation functions are 
measured at a rate depending on the speed of aging 
of different Laponite suspensions. The aging time for 
which the time-averaged correlation functions are not 
equal to their ensemble-averaged values, i.e. their val- 
ues change from one position to another in the sam- 
ple, defines the ergodicity-breaking point teb- For aging 
times ta > teb, we calculate the ensemble-averaged elec- 
tric field correlation function i. e. intermediate scattering 
function f{q,t) from the time-averaged intensity corre- 
lation function gtiq,t) and ensemble-averaged intensity 
If; measured by rotating the sample at different heights 

fiq,t) = 1 + iIt/lE){[gt{q,t)-gtiq,0) + - 1} (1) 

Figure ([1]) shows the evolution of intermediate scatter- 
ing functions f{q,t) for two different samples. In both 
cases correlation functions evolve from an ergodic state 
to a non-ergodic state (the correlation function no longer 
decays to zero), as the system ages. However, the low- 
and high-concentration samples are seen to behave in a 
distinctly different manner. In the non-ergodic regime, 
the aging rate of the system can be quantified by mea- 
suring the time evolution of the non-ergodicity parameter 
/((?, do), defined as the long-time limit of the intermedi- 
ate scattering function. Plotting the non-ergodicity pa- 
rameter as a function of scaled aging time reveals that its 
evolution in a range of samples with different Laponite 
concentrations collapses onto two distinct master curves 
when plotted (Fig. 2a) as a function of reduced aging 
time {ta/teb — !)■ In the first group of samples (low con- 
centrations) the non-ergodicity parameter almost reaches 
unity: the colloidal particles are completely blocked, sug- 
gesting that they are rigidly held in place, as they would 
be in a gel-like structure. In the second group (high con- 
centrations), the non-ergodicity parameter evolves at a 
slower rate and goes to ~ 0.8 at late times, indicating 
that there is still some freedom for the particles to move: 
the hindrance of the particles is only sterical, as it would 
be in a glass. In fact measuring the non-ergodicity pa- 
rameter for a couple of glassy samples until ta/teb = 15, 
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FIG. 1: Evolution of intermediate scattering function for 
two Laponite concentrations, at scattering angle 90°. The 
symbols present the measured correlation functions at in- 
creasing waiting times (from left to right) that are {ta = 
1, 49, 71, 80, 85, 89, 91, 99, 112, 141 days) for C = 0.2 wt% and 
{ta = 7, 40, 54, 71, 86, 113, 260, 1356 min) for C = 3.5 wt%. In 
both panel, the lines show the fits with Aexp(— t/ri) + (1 — 
A)eM-{t/r2f)- 



we found that f{q, oo) did not exceed 0.835. 

Perhaps the most striking observation is that for the 
intermediate concentrations 1.3 < C < 2.3wt%, identical 
samples may evolve at very different rates, thus having a 
very different ergodicity-breaking point. In Fig. ^jp), we 
have plotted the ergodicity-breaking time of the ensemble 
of the samples we have measured as a function of con- 
centration. The samples fall in two separate groups, and 
the samples in the intermediate concentration region fall 
in either of the two. This suggests that the intermediate 
concentration samples have two options, either following 
the same trend as the samples of high concentration or 
behaving similarly to the samples of low concentrations. 

The difference between the two states becomes clear if 
we look at the rest of the data. Before ergodicity break- 
ing, as can be observed from the fits shown in Fig. ([T]), 
both the low-and the high concentration samples can be 
described by a sum of a single exponential and stretched 
exponential A exp ( - t/n ) -h ( 1 - A) exp( - (t /ra )'' ) 0, i, i] . 
Ti is related to the inverse of the short-time diffusion 
Ti = l/(Dsg^); as is shown in fig. [3l the short time 
diffusion coefficient is almost a constant for the samples 
of high concentration (C > 2.3 wt%); this corresponds 
to the 'rattling in the cage' motion reported earlier for 
colloidal glasses 18]. However, Ds decreases significantly 
with aging time for low concentrations (C < lAwt%) 
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FIG. 2: a) The evolution of the non-ergodicity parameter 
/(g, co) versus reduced aging time ta/^eb — 1 for different 
Laponite samples. The colloid concentrations are shown in 
the legend. The samples can be divided into two groups ac- 
cording to the evolution of non-ergodicity parameters. In both 
figures (and the following ones), the open symbols correspond 
to gels, the filled symbols to the glass, b) The ergodicity- 
breaking time as a function of concentration. 
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FIG. 3: a) The evolution of short-time translational diffusion 
normalized to its initial value {ta ~ 0) as a function of ta/teb 



which is again related to the incorporation of the col- 
loidal particles in a gel network 9| . In addition, the slow 
relaxation time T2, is found to grow exponentially with 
aging time for high concentrations, in agreement with 
earlier observations for the glassy state [sffsil. However, 
for the gel phase, the relaxation time increases faster than 
exponentially (Fig.H^); this is likely to be related to the 
formation of small clusters in the beginning that subse- 
quently aggregate to form a macroscopic structure (with 
a large relaxation time) as in diffusion- limited clus- 

ter aggregation (DLCA). 

The measured static structure factor (Fig. O indeed 




FIG. 4; a) Slow relaxation time T2 as a function of scaled 
aging time b) non-ergodicity parameter as a function of scaled 
aging time. These samples behave like a glass at the early 
stages of aging and a after t^/t^ij ~ 1.2 they evolve according 
to the gel line. The filled lines show the glass and dotted lines 
show the gel line, obtained by smoothed averaging over all 
the samples measured. 



provides two further pieces of evidence for structure for- 
mation. First, for all the low-concentration samples, the 
scattered intensity consistently increases with aging time, 
while for the high concentration ones it does not evolve 
much (Fig. [5^) . The increase of scattered light intensity 
is usually attributed to formation a network or clusters of 
particles . Il4lll6j . Second, although the range of wave 
vectors is rather limited, measurement of the static struc- 
ture factor S{q) for the different samples , a power law 
behavior for S{q) appears to be observed for the lowest 
concentration samples with an exponent 1.1 ± 0.2. This 
exponent is indeed characteristic of DLCA, in which less 
compact clusters are formed than in diffusion-limited ag- 
gregation [14, 16]. To the contrary, an almost flat struc- 
ture factor is found for the highest concentrations, in- 
dicating homogeneity (Fig. (03), very characteristic of a 
glass [7[. The noise at low q for these measurements is 
probably due to imperfections of the measurement cell 
that scatter light at small angles. 

Combining all these data, we identify the low- 
concentration samples as colloidal gels and the high- 
concentration ones as colloidal glasses. Intermediate con- 
centrations can be either gels or glasses at late times, with 
no way of telling beforehand how the sample is going to 
evolve. That the path towards these non-ergodic states 
is complicated follows from the observation that samples 
may 'hesitate' for a long time between the two states, and 
may evolve in one direction to end up in the other. In 
Fig. m some of the 'hesitating' samples are indicated in 
color. It clearly shows how a few samples in the interme- 
diate concentration region that behaved consistently like 
glassy samples before the ergodicity breaking point, end 
up as gels at late times. Perhaps even more surprisingly. 
Fig. [4Jd shows that even if one looks at a single observ- 
able such as the non-ergodicity parameter, a crossover 
between the two behaviors can be observed; this is most 
evident in the data for 1.9wt%. The data shown in the 
figure are mere examples; approximately 20% of the sam- 
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FIG. 5: a) Scatted intensity at scattering angle 90° as a func- 
tion reduced aging time. So as to focus on the effect of aging, 
we have normalized the intensity to its value at the beginning 
of aging, b) The scattered intensity as a function of dimen- 
sionless wave vector qR for several Laponite concentrations. 
To be able to compare the q-dependence of different samples, 
we have normalized the data to their value at highest mea- 
sured qR — 0.38. As we observe the higher the Laponite 
concentration, the more homogenous it becomes. 



pies in the intermediate concentration region behaved in 
an ambiguous way in the sense that they seemed to have 
a hard time to 'decide' whether they were glasses or gels. 

In conclusion, the nature of the non-ergodic state in 
Laponite suspensions has been the subiect of consid- 
erable controversy: both colloidal gel [a, and col- 
loidal glass formation 0, [ll[ have been invoked and were 
thought to be mutually exclusive 0, • We have shown 
here that gel and glassy states of Laponite both exist and 
are well-defined in the limit of low and high concentra- 
tions. In gels the main cause of aging is the building up 
of a network. Slowing down of particle motion will hap- 
pen when the particles are trapped in the network, which 
may be formed already for very low particle concentra- 
tions. On the other hand, the slowing down of motion 
in glasses happens due to topological constraints jl^, 
that hinder the motion of particles. Therefore they are 
usually formed in samples of higher concentrations. We 
here provided clear experimental evidence for how the 
distinction of gel and glassy states can be made experi- 
mentally. 

This allows for the observation that for intermediate 
concentrations, the transition to non-ergodicity can oc- 
cur in either direction (gel or glass), and may be ac- 



companied by 'hesitations' between the two directions. 
A qualitative explanation for this behavior is provided 
by the free energy landscape picture of slow dynamics. 
Here, one attributes the slow structural relaxation to 
the complex pathways that connect the configurational 
states on the multidimensional potential surface. For 
Laponite our data suggest that there are at least two 
metastable minima in the free energy corresponding to 
gel and glass states, and that different pathways to- 
wards these non-equilibrium states exist, providing the 
first evidence for the existence of such a complicated free- 
energy landscape. This phenomenon should be general 
for colloidal systems with attractive interactions between 
the particles; indeed the recent discovery of "attractive 
glasses" for spherical colloids [S] also suggests that gels 
and glasses are not necessarily clearly distinct states of 
matter, but rather metastable minima in an otherwise 
complicated free-energy landscape, resulting from both 
steric and attractive interactions, as suggested by some 
simulations fv?\ and experiments [l6l |. 
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